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The spurious interaction of quantum systems with their environment known as decoherence leads, 
as a function of time, to a decay of coherence of superposition states. Since the interactions be- 
tween system and environment are local, they can also cause a loss of spatial coherence: correlations 
between spatially distant parts of the system are lost and the equilibrium states can become lo- 
calized. This effect limits the distance over which quantum information can be transmitted, e.g., 
along a spin chain. We investigate this issue in a nuclear magnetic resonance quantum simulator, 
where it is possible to monitor the spreading of quantum information in a three-dimensional net- 
work: states that are initially localized on individual spins (qubits) spread under the influence of a 
suitable Hamiltonian apparently without limits. If we add a perturbation to this Hamiltonian, the 
spreading stops and the system reaches a limiting size, which becomes smaller as the strength of the 
perturbation increases. This limiting size appears to represent a dynamical equilibrium. We present 
a phenomenological model to describe these results. 

PACS numbers: 03.65.Yz, 03.67.Ac, 72.15.Rn, 76.60.-k 



I. INTRODUCTION 

Controlling quantum mechanical systems has received 
increasing attention in recent years pQ, mainly because 
we are starting to be able to perform computations with 
quantum mechanical systems. This has the potential of 
solving computational problems for which no efficient al- 
gorithm exists for classical computers [2H3]. However, we 
will only be able to realize this potential if we learn to 
control large quantum systems with high reliability. The 
control of small quantum systems has been thoroughly 
explored in the last years pQ , however the study and con- 
trol of large quantum systems has not been tackled so 
far. The reason for this is partly the difficulty of sim- 
ulating large quantum systems on classical computers, 
which limits the number of qubits to about 20 when the 
system is in a pure state [5j [6] . Additionally, it is also 
possible to calculate the dynamics of mixed states if the 
initial state is localized by using quantum parallelism of 
a single pure state evolution [8]. On the other side, 
the present primitive state of quantum computation only 
allows incomplete control of large quantum states. So 
far, the only physical system that offers this possibility 
is nuclear magnetic resonance (NMR) of dipolar coupled 
spins [9HTT]. In particular, the main problems are the 
lack of individual addressing of qubits and decoherence. 
The latter degrades the quantum information of a given 
state [12] . Decoherence was shown to increase as the 
size of the quantum system increases, making the largest 
systems the most sensitive to perturbations [9HTTj [T~3l 
[T5] . While reducing the perturbation strength typically 
reduces decoherence, in complex and large systems, the 
decoherence time can be independent of the perturbation 
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strength over a range of coupling strengths p~6j H7] . 

Decoherence is well known as a process that causes the 
decay of quantum information. Avoiding or reducing de- 
coherence is thus the main ingredient for implementing 
large scale quantum computers. Several techniques have 
been proposed for this purpose, including dynamical de- 
coupling [18], decoherence- free subspaces [19], and quan- 
tum error correction [20 , 21 j. These proposals have been 
tested on small systems of nuclear spins [22] , trapped ions 
[23] or spin model quantum memories [24] [25] . Some of 
these techniques were successfully applied to large quan- 
tum systems with thousands of qubits where their de- 
coherence time was extended by almost two orders of 
magnitude [10] [15] . 

Decoherence not only affects the survival time of quan- 
tum information, it also affects the distance over which it 
can be transmitted [26-34 . For example, spin systems, 
and in particular spin chains, can be used to transfer 
quantum information over large distances [35, 36 . This 
kind of systems were studied with liquid state NMR for 
small numbers of spins [37H4T] and slightly larger num- 
bers by solid-state NMR [42H44] . However, once deco- 
herence is considered, in the simplest two spin quantum 
channel, it was shown that when the effective system- 
environment interaction exceeds a given strength, it be- 
comes impossible to perform even a simple SWAP oper- 
ation [33] [45] . Instead of having the expected oscillatory 
transfer of a state going forth and back between the two 
spins, an overdamped dynamics due to the appearance of 
a localized state appears at a critical value or exceptional 
point of the perturbation strength [45 -49 . Similarly, this 
could be observed in a 3-spin chain. If one spin is suffer- 
ing a perturbation, when it exceeds a given strength the 
dynamics localizes in the remaining two spins [50 . In a 
more general situation with longer spin-chains, it was re- 
cently pointed out that imperfections or disorder of the 
spin-coupling that drives the state transfer can induce 
localization of the quantum information [26 - 31, 34 in a 
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process related to Anderson localization [5TJ [52] . In more 
complex 3D spin-network topologies, we demonstrated 
experimentally a similar behavior by studying the local- 
ization effects induced by the finite precision of quantum 
gate operations used for transferring quantum states [32] . 

In this paper, we extend our previous work [32], where 
we prepared the system first as individual, uncorrelated 
spins and measured the build-up of clusters of correlated 
spins of increasing size. Introducing a perturbation to the 
Hamiltonian that generates these clusters, we find that 
the size of the clusters reaches an upper bound. This up- 
per bound appears to be a dynamic equilibrium: if the 
cluster size is initially larger than this equilibrium value, 
it decreases under the effect of the perturbed Hamilto- 
nian, while the unperturbed Hamiltonian leads to an in- 
crease. The equilibrium size decreases with increasing 
strength of the perturbation. For these experiments, we 
use a solid-state NMR quantum simulator. While this 
system does not allow addressing of individual qubits, 
it represents an excellent test-bed for studying different 
aspects of decoherence and information transfer. 

II. THE QUANTUM SIMULATOR 



We performed all experiments on a home-built solid 
state NMR spectrometer with a X H resonance frequency 
of uj z /2tt = 300 MHz. The spins are the protons of 
polycrystalline adamantane, where the average strength 
of the average dipolar interaction, determined from the 
width of the resonance line, is 7.9 kHz. 

III. GROWTH OF SPIN CLUSTERS 

A. Cluster size 

The spin clusters that we consider can be written in 
product operator form as 

Here, the indexes Z, o,p identify the spins involved in the 
given cluster. We write K for the number of terms in 
this product, i.e., for the number of spins in the cluster. 
Experimentally, we generate these clusters using an NMR 
method developed by Pines and coworkers [54] [55] . It is 
based on generating an average Hamiltonian j-Lo that does 
not commute with the thermal equilibrium state 



We consider a system of equivalent spins I = 1/2 in 
the presence of a strong magnetic field. The Hamiltonian 
of the system is 



i i<j 



2iiii-(r x p x + r y p y ) , (i) 



where uj z is the Larmor frequency, dij the coupling con- 
stants and I l x ,Iy, l\ are the spin operators that can be 
represented by P u = \o u with a u Pauli operators. 1~L Z 

represents the Zeeman interaction and 1-Ldd the dipolar 
interaction [53] . The latter is truncated to commute with 
the strong Zeeman interaction (uj z ^> d^-), assuming that 
the effects of its non-commuting part are negligible. In a 
frame of reference rotating at the Larmor frequency [53] , 
the Hamiltonian of the spin system reduces to l~Ldd- 

The quantum simulations start from the high- 
temperature thermal equilibrium [53 . Using the nota- 
tion I z = ^ I\ , we can write the thermal equilibrium 
state as 



Po 



1 



huj z 



h /Tr{i}. 



(2) 



The unity operator 1 commutes with all operators, in- 
cluding the Hamiltonian and the density operator. It 
does not contribute to the observable signal and it is 
therefore convenient to exclude it from the density oper- 
ator. The resulting density operator for the initial state 
of the system is then po oc I z . In this state, the spins are 
uncorrelated and the density operator commutes with the 
Hamiltonian 1-Ldd- 



u = -J2dij [rJi-im 



i<j 

i<j 



P+P+ + p_Il 



(3) 



where I± 



^/■jdzif-j^. In the usual computational 
basis, we write the states as |M, um) where M is the 
total magnetic quantum number, i.e., I z \M,tim) = 
M |M, tim) ? an d Um distinguishes different states with 
the same M. The Hamiltonian ([3| flips simultaneously 
two spins with the same orientation. Starting from the 
thermal equilibrium state, it generates a density operator 
where only elements p^ with 

AM = M(i) - M(j) = 2n, n = 0, 1, 2 . . . 

are populated. The index i and j refer to the different 
spin states in the Zeeman basis as described above. Such 
a density operator element pij is called a AM quantum 
coherence. Off-diagonal elements with AM = repre- 
sent zero-quantum coherences and diagonal elements cor- 
respond to populations. 

Figure [l] shows a graphical representation of the den- 
sity matrix in the computational basis. The black diago- 
nal line represents the population elements of the density 
matrix. The diagonal blocks (gray) represent coherences 
of p with AM = and moving away from these diagonal 
blocks, the AM multiple-quantum coherence (MQC) or- 
der increases and they are represented by different gray 
values. The solid arrows indicate the effect of 1~Lo on 
the thermal equilibrium density operator. Adding the el- 
ements of p with a given AM, we obtain a distribution 
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Figure 1: Density matrix in the computational basis. The 
black diagonal line represents the populations (Pop. in the 
legend). The diagonal blocks (dark gray) represent coherences 
of p with AM = (ZQC) and the non-diagonal blocks are the 
AM multiple- quantum coherence (MQC) blocks (gray tones). 
The solid and dashed arrows indicate the effect of Tio and the 
perturbation interaction E = Tidd respectively on the thermal 
equilibrium density operator. A schematic representation of 
a typical distribution is given by the Gaussian-like shape on 
the top-left corner. 



function of the MQC elements of the density matrix pam- 
A schematic representation of its distribution is given by 
the Gaussian-like shape on the upper-left corner of Fig. 
[T] This distribution is initially a delta function on the 
diagonal of p and spreads with time. 

To determine the average number of correlated spins, 
the technique relies on the fact that in a system of K 
spins, the number of transitions with a given AM follows 
a binomial distribution [56j [57] 



n (AM, K) = 



(2K)\ 



(K + AM)! (if —AM)! ' 



(4) 



For K ^> 1, the binomial distribution can be well approx- 
imated with a Gaussian 



n (AM, K) oc exp 



AM 2 
K 



(5) 



whose half-width at e is a = yK . Thus, we can 
determine the effective size of the spin clusters in a given 
state by measuring the distribution of the MQCs of its 
density operator p as a function of the coherence order 
AM. 



B. Growth 

Figure [2] shows two distribution functions of the MQCs 
for two different durations Ntq of the evolution under 
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Figure 2: (Color online) Multiple-quantum coherence (MQC) 
distributions for two different evolution times. The cluster- 
size is extracted from the variance a = \[K. 



Ho. They clearly demonstrate the increasing width of 
the MQC distribution and thereby the increasing size of 
the spin clusters. 

The elements of the density operator with a given AM 
can be distinguished by rotating the system around the 
z— axis. A rotation <\> z = e~ l ^ lz by <p changes the density 
operator as 



p (</>, t) = kp (t) 4>\ = £ feM (*) e ' AM *. 

AM 



(6) 



where p(t) = e~ %Hot poe %Hot and /?am(£) contains the ele- 
ments that involve coherences of order AM. The experi- 
mental observables in NMR are I x , I y , and by means of a 
7r/2 pulse, also I z . They are operators with AM = ±1 or 
AM = 0, respectively. Thus we are not able to measure 
all the elements of the density matrix. In order to quan- 
tify the p AM blocks, we use the pulse sequence shown in 
Fig. H 

The system initially evolves for a period of duration 
Nto under the Hamiltonian (j~Lo^ = fizHofit, ^.e., 

Po > Pcf> (Nr ) = (f) z p (Nr ) <j>\ 

= z e-^ Nr °p o e™° Nr °ft z 



AM 

£/W^r )e^. (7) 
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Figure 3: (Color online) NMR sequence for generating large 
spin clusters. The effective Hamiltonian Tio is generated by 
the sequence of ty/2 pulses shown in the upper part of the 
figure, where A' = 2 A + r p and r p is the tt/2 pulse duration 
[55]. 
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A subsequent evolution period with the same duration 
Ntq under —Ho, which causes an evolution backward in 
time, yields the final density operator 



(-Ho)Nto 



> Pf (2Nr ) 



JHqNto 



Pcj, (Ntq) < 



-iUoNro 



= J2[e iiloNTo PAM(Nr )e- i ^Nr 

AM 



iAMcf) 



■ (8) 



If I z is the NMR observable, we obtain the signal 
S( ( p,NTQ)=Tv{i z p f (2Nt )} 

= Tr {e- in ° NT0 poe™ ^ fa (iVr )} 

= Tr{p(NT )fo(NT )} 

= J2e^ M Tr{pl M (Nr )} (9) 



AM 



^e^ AM A(AM) ; 



(10) 



AM 



where A(AM) are the amplitudes of the MQ spectrum. 
To extract these amplitudes from the experimental data, 
we measure the signal S (0, Ntq) as a function of (j) and 
perform a Fourier transform with respect to <fi. We 
quantify the cluster size K by inverting the relation 
A (AM) = n (AM, K) (see Ref. [55] for details). 

If the system evolves under the Hamiltonian (J3|, the 
width of the A (AM) distribution increases indefinitely, 
as shown in Fig. [4j The main panel shows the MQC 
distributions A (AM) for different evolution times. We 
quantify the spreading of the MQC distribution by mea- 
suring their widths (a) at A (a, 7Vro)/A(0, Ntq) = 1/e for 
different evolution times. The inset of the figure shows 
the time evolution of a. Equivalently, its corresponding 
cluster size of correlated spins also grows indefinitely, as 
shown in Figure [5] 

The Hamiltonian of Eq. ^ is prepared by means of 
the standard NMR sequence p4[ [55] shown in the upper 
part of Fig. [3j This sequence consists of 7r/2 rotations 



Figure 4: (Color online) Time evolution of the multiple quan- 
tum coherence distribution. The main panel shows the time 
evolution of the MQC spectrum A(AM). The inset shows the 
time evolution of a, where A(a, Ntq)/A(0, Ntq) = 1/e. 
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Figure 5: (Color online) Time evolution of the cluster size 
of correlated spins with the Hamiltonian %q (black square 
points) . The main panel shows it in log scale manifesting the 
region where the growth is exponential and the inset shows 
it in log- log scale manifesting its power law growth for long 
times. 



of the spins separated by periods of free precession, A 
and A 7 , under the dipolar Hamiltonian Hdd an d it was 
shown to approximate quite well the ideal Hamiltonian 
7-Lq. It is inverted experimentally by shifting the phase 
of all RF pulses by ±7r/2 [54 . We used pulse durations 
t tt/2 = 2.8/is and delays A = 2/is and A' = 2A + t p 
giving a cycle time tq = 57.6/is. 
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C. MQC distributions and cluster-size 

In this section we summarize previous works that de- 
scribed the MQC distribution, its evolution and how it is 
related with the cluster-size of correlated spins in order to 
compare their results with our experimental observation. 
So far, it has not been possible to derive a consistent 
theoretical model of the processes involved in this type 
of many-body system. The most accepted and simplest 
model for interpreting the MQC distributions was pro- 
posed by Baum et al. [55]. It assumes equal excitation 
probabilities (for a given system size) for every coherence 
term of the density matrix and therefore predicts a Gaus- 
sian distribution whose variance is related to the number 
of correlated spins. This is the usually adopted model 
because it provides the simplest qualitative description 
of the MQC evolution. However, the experimentally ob- 
served distributions are not always Gaussian [58]. This 
is also the case in our experiments: the data shown in 
Fig. [4] are better described by an exponential distribution 
A(AM). Our results are therefore more consistent with 
the distribution discussed by Lacelle et al. [58 , who pre- 
dicted an exponential MQC distribution and found agree- 
ment with earlier experimental results from adamant ane. 
Additionally, they predicted a dynamical scale invariance 
in the growth process of the MQC. 

Due to the difficulty of treating many-spin systems rig- 
orously, quantum mechanical predictions for the MQC 
profile shapes and their interpretations are still missing. 
For example, if the MQC profile is not a Gaussian distri- 
bution, it is not trivial how to extract from it the cluster- 
size. While, most of previous works studying the distri- 
bution of MQCs have assumed Gaussian distributions, it 
is more convenient to use a parameter that is indepen- 
dent of the MQC distribution form. Khitrin proposed to 
use the second moment of the MQC spectrum for this 
purpose [59 . This parameter was shown to be related 
to the second moment of the system Hamiltonian, in our 
case given by Eq. ([3| and K(t) is proportional to the 
second moment of the MQC distribution. 

To determine the cluster-size from our experimen- 
tal data we measured the half-width a = \/K of the 
MQC distributions at 1/e. We chose a = \/K be- 
cause an exponential MQC distribution of the form 

exp ^— AM/ ^K(t)^j has the same width as the Gaus- 
sian distribution ([5|. Assuming this distribution, its sec- 
ond moment is 2K while for a Gaussian distribution it 
is K/2. With the exception of the Gaussian model for a 
MQC distribution, there is no rigorous model that pro- 
vides the exact factor that converts the second moment 
of the distribution to the cluster size K. However, they 
are of the same order of magnitude and thus the error of 
the cluster-size determination is a scale factor of order 1. 
In our experiments, the MQC distribution has the same 
shape during the complete range of the experimentally 
accessible evolution time. Accordingly, any error of the 
determination of the cluster-size is thus a constant factor 



independent of the value of K and this does not change 
the conclusions of this article. 

The growth depends on the spin-coupling network 
topology as was observed experimentally [55j I58H63] . 
Mainly three types of growth were observed: (i) indefinite 
growth, (ii) localized growth and (iii) localized growth at 
a first stage (intra-molecular-like localization) and then 
a further indefinite growth (inter- molecular). 

Results for adamantane fall into the first category. 
Some models suggest a power law for the growth [61] [64] 
and these preliminary results [58 , 61 seem to match with 
an effective 3D spin-coupling network topology. That 3D 
behavior seems to be achieved only for cluster sizes above 
1000 spins [6T] , However, more recent works suggest an 
exponential growth of the cluster size when many neigh- 
boring spins are contained in the spin-coupling topology 
[62] [63] . This is the case of adamantane and the experi- 
ments agree well with that prediction. 

Experimentally, the data shown in Fig. [5] suggest three 
stages of the cluster size evolution: (i) An initial period 
where the evolution cannot be described as exponential 
or power law, (ii) a period of exponential growth, and (iii) 
a power law regime. During the first stage, the MQC dis- 
tribution changes from a Gaussian-like distribution to an 
exponential distribution. Thus, the determination of the 
cluster size could be affected by two kinds of systematic 
errors: (i) at small cluster sizes, the small number of cor- 
related spins may make the statistical assumptions ques- 
tionable and (ii) the change of distribution from Gaus- 
sian to exponential may change the appropriate scale fac- 
tors. The exponential growth agrees with the behavior 
expected by theoretical predictions for this kind of sys- 
tem [62] [63] . And finally the power law behavior could 
be due to the following two reasons: (i) as predicted by 
Lacelle [61] after a certain number of spins (around 1000 
in adamantane) the effective spin-network topology turns 
into a 3D spin-coupling network, which leads to a power 
law growth or (ii) because the experimental generation 
of Ho contains non-idealities, they can be considered as 
a perturbation that generates localization effects, as was 
demonstrated in Ref. [32 and discussed in the following 
sections. Something similar to the latter point could be 
interpreted from the predictions of Ref. [63] . 



IV. EFFECT OF PERTURBATION 

A. Perturbed evolution 

The evolution under the Hamiltonian Ho can be re- 
versed completely by changing the Hamiltonian from Ho 
to —Ho- If this time reversal is perfect, the signal ([9| for 
phase (j) = is independent of the evolution time iWo, 
^2 AM A(AM,Nr ) = const. 

Experimentally, the Hamiltonian Ho is generated as 
an effective Hamiltonian. Because of experimental im- 
perfections, it always deviates from the ideal Hamilto- 
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Figure 6: (Color online) NMR sequence for the quantum sim- 
ulations. A perturbed evolution is achieved when te ^ 0. 
The effective Hamiltonian Ho is generated by the sequence of 
7r/2 pulses shown in the upper part of Fig. Inland XI = Hdd is 
the free evolution Hamiltonian. 



Taking into account now the complete sequence of evo- 
lutions given by Fig. [6] 

(n eii )Nr c A 

Po Jj -^ pf* (Nt c ) = fopF* (Nr c ) ^ 



AM 

Ep2S(JVt c )c^*(13) 



AM 



nian (3]). As a result, the actual dynamics deviates from 
the ideal one and, in particular, we cannot exactly in- 
vert the perturbed Hamiltonian and thus revert the time 
evolution. Because of the deviations in Hq : the quantity 
Sam A(AM) is no longer conserved, but decays with in- 
creasing evolution time. This decay is not uniform, but 
it affects mostly those components of the density oper- 
ator that correspond to strongly delocalized coherence. 
As a result, the spreading of the information is attenu- 
ated and the system becomes effectively localized. It is 
this latter effect that we want to study here; we isolate 
it from the overall decrease of the signal by normalizing 
the MQ spectra such that the total signal Eam A(&M) 
for (j) = is again constant in time. 

To analyze this deviation from the ideal evolution, we 
introduce a perturbation E, whose strength we can con- 
trol experimentally and study the behavior of the system 
as a function of the perturbation strength. We choose the 
dipole-dipole coupling for this perturbation, E = Hdd-, 
which is a local interaction: every spin interacts mostly 
with its nearest neighbors, while the coupling strength 
with more distant spins drops off as 1/d 3 . 

We add this Hamiltonian to the ideal Hamiltonian Ho 
by concatenating short evolution periods under Hdd with 
evolution periods under Ho- We label the durations of 
the two time periods and To, as shown in Fig. [6j 

When the duration r c = ro + 7e of each cycle is short 
compared to the inverse of the dipolar couplings d^-, 
the resulting evolution can be described by the effective 
Hamiltonian 



?4ff = (1 -p)W>o +pE, 



(11) 



where the relative strength p = /r c of the perturbation 
£ = ^dd can be controlled by adjusting the duration T£. 
In the experiment, we compare the artificially perturbed 
evolution of H e ff with the Ho evolution with its intrinsic 
errors. While the intrinsic errors reduce the signal or the 
overall fidelity, they do not cause localization on the time 
scale of our experiments (see Fig. [5]). 

Considering now this perturbation, starting from ther- 
mal equilibrium, the state of the system at the end of TV 
cycles is 



(-H )Nt nj , 

± } - >pf* (Nt c + Nt ) = 

= e iHoNr -n ei{ ( Nr j e -in Nr 

= [e tn ° Nro PAM (Nt c ) e~^ N ^ 

AM 



JAMcf) 



■ (14) 



(Nt c ) 



-iV. e ffNr c ^ e in e ffNr c 



(12) 



Thus the NMR signal, which is measured after the last 
backward evolution exp jzHoA/Yo j, can be written as 

S (</>, Nt + Nt c ) = Tr {i z p** (Nt c + Nt ) } 

= Tr je-^oWro^giWoA^Heff (JVr c )| 

= Tr{p«» (Nt ) pf" (Nt c )} 

= E e ^ AMTr {pa°m (Nto) p^t! (Nr c )} 

AM 

(15) 

= ^e^ AM i(AM). (16) 

AM 

Reinterpreting these expressions, we define the effec- 
tive observable A = e ~ i1loNTo p e i1loNTo = p Ho (Ntq), 
where p n ° is the density operator of the unperturbed 
evolution. The NMR signal becomes then S((j), Nr c ) = 

Tr{Apf«(Nr c )}. 

For the ideal evolution (p = 0), Eq. ^ is recovered, 
where A(AM) correspond to the squared amplitudes of 
the density operator elements p^° M (Ntq) with coherence 
order AM. For the perturbed evolution, [p ^ 0), they 
are reduced by the overlap of the actual density opera- 
tor elements /^m i^ r c) with the ideal ones. We extract 
these amplitudes by performing a Fourier transformation 
with respect to (j). Figure [7] shows, as an example, the re- 
sulting A(AM) as a function of time for p = 0.108. The 
main panel represents the MQC distributions for different 
evolution times. The main difference compared to Fig. [4] 
is that the MQC spectrum does not spread indefinitely 
[32] . We consider this saturation of the spreading of the 
MQC spectrum evolution as evidence of localization due 
to the perturbation (see below). The localization effects 
are easily visualized by directly comparing the generation 
of high-order multiple quantum coherences of the unper- 
turbed case of Fig. [4] with the perturbed one of Fig. [7[ 
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Figure 7: (Color online) Time evolution of the multiple quan- 
tum coherence distribution under a perturbation p = 0.108. 
The main panel shows the time evolution of the MQC spec- 
trum A(AM). The inset shows the time evolution of the stan- 
dard deviation a for the unperturbed case compared with the 
perturbed case. This comparison shows directly the localiza- 
tion. 



While the distribution spreads continuously in Fig. [4j it 
reaches a limiting value in Fig. [7[ The inset of Fig. [7] 
shows the time evolution of the width a for p = 0.108 
compared with the unperturbed case. 

From the width of the MQC spectrum, we calculate 
the size of the spin clusters. The black squares of Fig. [8] 
show the average number of correlated spins as a func- 
tion of time for an unperturbed evolution, p = 0. This 
is the log scale representation of the curve of Fig. |5j 
The other symbols of the figure show the evolution of 
the number of correlated spins for different values of p. 
Initially, the cluster size K(Nr c ) starts to grow as in the 
unperturbed evolution, but then it saturates after a time 
that decreases with increasing perturbation strength p. 
As we explained above, we consider this as evidence of 
localization induced by the perturbation and it is related 
to spatial localization. The size of the cluster at which 
this saturation occurs is also determined by the strength 
of the perturbation: increasing perturbation strength re- 
duces the limiting cluster size. 



B. Interpretation 



The terms in Eq. (15) are related to the fidelity of 
the corresponding density operator component pam with 
respect to the corresponding component resulting from 
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Figure 8: (Color online) Time evolution of the cluster size for 
different perturbation strengths. The black squares represent 
the unperturbed time evolution and the other symbols corre- 
spond to the perturbed evolutions according to the legend. 



the unperturbed evolution 



/am = 



^{p!°m (Nr )f%° M (Nt )} 



(17) 



which reaches unity for vanishing perturbation (p = 0). 
We therefore consider the reduction as a quantitative 
measure of the effect of the perturbation. The decoher- 
ence of different AM-blocks of the density matrix was 
studied when a given quantum state prepared by evolu- 
tion under Ho evolves under a pure Hdd — S perturba- 
tion [H[T3]. However, the quantification and characteri- 
zation of how a perturbation disturbs the different blocks 
of the density matrix during the creation of spin-cluster 
states is still not known. This merits further studies like 
in Ref. [65] . Here, we focus only on the effect on the num- 
ber of correlated spins and do not consider the change of 
the overall amplitude. Thus we normalize the integral of 
each spectrum and determine its width by fitting it to a 
Gaussian or an exponential. 

The spreading of the MQC spectrum is generated by 
the effective Hamiltonian created with the sequence of 
Fig. [l] The cluster size that we determine here corre- 
sponds to an overlap of the actual state with the ideal 
state resulting from unperturbed evolution. This is simi- 
lar to the fidelity measure in quantum computing, where 
the agreement between the actual state of the system 
with the target state is measured. In our case, the tar- 
get state is a growing cluster, while the actual state grows 
only for some time until it reaches a limiting size - typical 
for localization. The size of this localized state decreases 
with increasing perturbation strength. 

During our experiment, the magnetization is uniform 
throughout the sample, so the process does not lead to 
a spatial redistribution of magnetization. However, since 
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Figure 9: (Color online) NMR pulse sequence for preparing 
different initial clusters sizes and subsequently evolving them 
in the presence of a perturbation. 



we measure the number of correlated spins, we can at- 
tribute a length scale to the resulting state. Given a 
suitable initial state, the same process would generate 
highly entangled multi-qubit states. In the presence of a 
perturbation, the number of qubits that can be entangled 
in this way is limited to the size of the resulting cluster. 



V. EVIDENCE FOR A DYNAMICAL 
EQUILIBRIUM CLUSTER SIZE 

The experimental results presented above show that 
the cluster size reaches a stationary value. It remains to 
be seen if this limiting size results from a slow-down in 
the growth [29] or it represents a dynamic equilibrium 
state. We showed that actually they achieve a dynamical 
equilibrium state [32] . In order to do that we repeated the 
previous experiment for a series of initial conditions cor- 
responding to different clusters sizes. Figure [9] shows the 
corresponding pulse sequence: The initial state prepara- 
tion, consisting of an evolution of duration TVoTo under 
the unperturbed Hamiltonian T-Lq, generates clusters of 
size Kq. During the subsequent perturbed evolution of 



duration 7Vr c , these initial clusters evolve and Eq. (15) 
becomes 

S (</>, Nr c ) = e^ AM A(AM) 

AM 

E e #AM Tr {p% (Nt , N t ) fftfi (Nr c , N t )} . 

(18) 



AM 



Figure [10| shows the results from these measurements. 
The symbols represent the amplitudes A (AM) of the dif- 
ferent multiple quantum coherences as a function of the 
coherence order AM for different evolution times. Start- 
ing from an initial cluster size Kq = 141 for two different 
perturbations strengths (p = 0.080 in panel a, p = 0.025 
in panel b). In this figure, the A(AM) values are not 
normalized. We find that the width of the MQC distri- 
butions contracts as a function of time in panel a), but 
expands in panel b). Like in the case where we start from 
size Kq = 1, this expansion does not continue indefinitely, 
but it saturates. Similarly, the contraction also reaches 
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Figure 10: (Color online) Time evolution of the multiple quan- 
tum coherence distribution in the presence of a perturbation 
starting from an initial sate with Ko c± 141 correlated spins. 
The initial MQC distribution of the initial state is given by the 
black squares and other colored symbols gives its further time 
evolution for two different perturbation strengths p = 0.080 
(Panel a) and p = 0.025 (Panel b). 



an equilibrium value. This is evident at the longest evolu- 
tion times, where the gradients of the MQC distributions, 
which give the cluster sizes, are parallel in the semi-log 
scale representation. 

To quantify this, we determined the width of the MQC 



spectra and from that the cluster-size. Figure 11 shows 
the evolution of the cluster-size when we start from 
the initial value Kq = 141 for different perturbation 
strengths (empty symbols). The solid symbols represent 
the cluster-size evolution starting from Kq = 1, i.e., the 
curves shown in Fig. [8] For all perturbation strengths, 
the time series for a given perturbation strength converge 
to the same limiting value, independent of the initial con- 
dition. Thus, we find that the evolution leads to a limit- 
ing cluster size that varies with the perturbation strength 
p, but does not depend on the initial condition Kq. The 
limiting cluster sizes are represented by the symbols on 
the K — p plane in Fig. [TT] 

This is again confirmed by the results shown in Fig. 
[12} which summarizes the results for additional initial 
values and two perturbation strengths, p = 0.034 and 
p = 0.065. The filled symbols correspond to uncorrelated 
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Figure 11: (Color online) Time evolution of the cluster-size 
starting from two different initial states as a function of the 
perturbation strength. Filled symbols are evolutions from an 
uncorrelated initial state (Kq = 1) and the empty symbols 
start from an initial state with Kq — 141 correlated spins 
(marked with black squares joined with a solid line). The 
equilibrium cluster sizes are represented by symbols in the 
K — p plane. 



initial states (Kq — 1) and the empty symbols to various 
initial cluster sizes. For a given perturbation strength, 
the size of the spin clusters tends towards the same lim- 
iting value, independent of the initial condition. 



VI. LOCALIZATION SIZE VS. PERTURBATION 
STRENGTH 

A. Experimental evidence 

According to Figs. |8j [TT] and [l2j the size of the dy- 
namical equilibrium clusters decreases with increasing 
strength of the perturbation. For a quantitative anal- 
ysis of this dependence, we determined the size of the 
localized clusters from the data shown in Figs. |8} pT]and 
12 and plotted them against the perturbation strength in 
Fig. [l3j The diagonal line in the figure represents a linear 
fit to the experimental data (yellow/light line) that gives 
-Kioc °° p-i-86±o.05^ j e ^ ^ e c i us ter size decreases almost 
proportionally to the square of the perturbation strength. 
The error of the fit is indicated by the line-width. The 
limiting value for p = 1, K\ oc « 1, is consistent with 
the expectation that the system becomes completely lo- 
calized if the perturbation strength is significantly larger 
than the unperturbed Hamiltonian. 

The figure also summarizes the evolution of the cluster 
size before the static (localized) size is reached: If the ini- 
tial size is larger than the stationary value for the given 
perturbation strength, Kq > K\ oc , the cluster shrinks [in- 



Figure 12: (Color online) Time evolution of the cluster size 
starting from different initial sates. Filled symbols are evolu- 
tions from an uncorrelated initial state for two different per- 
turbation strengths given in the legend. Empty symbols start 
from an initial state with Kq correlated spins. 
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Figure 13: (Color online) Localized cluster size K\ oc (square 
symbols) of correlated spins versus the perturbation strength 
p. Three dynamical regimes for the evolution of the clus- 
ter size are identified depending of the number of correlated 
spins compared with the perturbation dependent localization 
value: a) a cluster size decreases, b) a cluster size increases c) 
stationary regime. 



set (a) in the figure, above the diagonal]. If it is smaller, 
Kq < K\ oc1 the size increases with time [inset (b), below 
the diagonal]. 



B. Phenomenological model 

A theoretical model that describes this behavior would 
be highly desirable, but is beyond the scope of this pa- 
per. Instead, we describe here a simple phenomenological 
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model that summarizes the observed behavior. 

For this purpose, we write the amplitudes of the MQC 
spectra as 



(19) 



The exponential form agrees with the experimental data 
(see, e.g. Figs. [4] and 10). According to Fig. |FJ the 
growth of the cluster size is exponential during a large 
part of the experiment. It may then be described by a 
differential equation of the type 



dK_ 
~dt 



(20) 



where a is the growth rate, which is proportional to the 
second moment of the NMR resonance line [62] . 

To derive an equation of motion for the perturbed 
evolution, we first consider the effect of the perturba- 
tion alone. As we have shown before [9], the interaction 
with the environment causes a decay of the MQC am- 
plitudes that depends on the MQC order AM. It was 
shown in Ref. [9] that over the range of our measure- 
ments the decay rate is almost ex AM. The decay of 
the MQ amplitudes during a time r is thus SA(AM) oc 
— A(AM)AMpbr, where p is the perturbation strength 
introduced above and b is the decay rate of single quan- 
tum coherences under the full perturbation, p = 1. 

During a short time r, the MQC spectrum evolves thus 
from (19) to 



A(AM,t + r) = Ae 



-AM 



(l/y/K^t+^+pbr) 



Here, K\ = Kq(1 + ar) is the size that the cluster size 
would reach, starting from K$ in the absence of the per- 
turbation. We can rewrite this as 



with 



A(AM, t + r)= 4e- AM /v^ 



K (l + ar) 



K (1+cxt) 



1 + 2y/K (l+aT)pbT + K (l + ar)(pbT) 2 
where for short r 



, (21) 



K (l + ar) 



l + 2^K (l + ar)pbT 

K (l + ar)(l - 2y/K (l + aT)pbr). (22) 



We now look for the stationary solution where K[ = Kq : 
K[=K = K (l + ar)(l - 2^K (1 + ar)pbr). (23) 

If we consider only terms of 0(r) in the infinitesimal time 
r, this is equivalent to 

(a - 2^K~vpb) = 0. (24) 
Solving for Kq, we find the stationary cluster size as 



a 

2bp 



(25) 



Considering the simplicity of this phenomenological 
model, this result agrees reasonably well with the ex- 
perimentally observed behavior. 



VII. DISCUSSION AND CONCLUSIONS 

As a step towards improved understanding of the evo- 
lution of large quantum systems, we have studied the 
spreading of information in a system of nuclear spins. 
Starting with single qubits, the information can spread 
to clusters of several thousand qubits. While decoherence 
is well known to limit the time for which quantum infor- 
mation can be used, we focused here on its effect on the 
distance over which a quantum state can be transferred. 
For that purpose, a locally stored state was left to evolve 
in a 3D spin-coupling network. Using standard NMR 
techniques we generated a Hamiltonian Ho that spreads 
the information and allows one to quantify the size of 
the resulting cluster of correlated spins. By comparing 
this unperturbed evolution with the evolution where a 
perturbation Hamiltonian of variable strength is added 
to Ho , we showed that the information becomes local- 
ized on a distance scale that decreases with increasing 
perturbation strength. 

Our experimental results demonstrate that a common 
dynamic equilibrium size of the localized state is achieved 
independent of the initial state consisting of different 
numbers of correlated spins. We developed a phenomeno- 
logical model to describe these effect, which may be at- 
tributed to the competition between the spreading evo- 
lution driven by Ho and its systematic reduction due to 
decoherence induced by the perturbation. These results 
are related with theoretical predictions on similar sys- 
tems that indicate a slow down of the spreading [29] and 
can induce localized states with a finite cluster-size [66] . 
Further extensions of these previous works could show if 
a dynamical equilibrium independent of the initial state 
appears as in our experiments. These previous works 
have considered disorder as the perturbation that pro- 
duces Anderson localization. 

The results presented here showed a transition in the 
spin dynamics from an indefinite spreading to a localized 
dynamics. A spin far away from the spin where the initial 
condition is stored would receive excitation at some time 
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if the perturbation is below a critical value, however, if 
the perturbation exceeds this threshold no excitation will 
arrive at this site. Our experiments show a transition for 
the cluster-size dynamics: if the cluster-size of the initial 
state exceeds the localization value, the it shrinks until it 
reaches the equilibrium value, but when the initial cluster 
size is lower than the localization value, it grows. 

These results may also be connected to our earlier find- 
ings that the decoherence rate of quantum states with 
many correlated qubits increases with the size of the sys- 
tem [9] as \[K, indicating that larger systems are more 
sensitive to perturbations. This increasing of the deco- 



herence rate, as the system size increases, balances the 
tendency of the system to spread. 
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